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I. INTRODUCTION 


Ave FECKGRCUNDT AND PURPOSE 

Curing recent years,there have been a number of scien- 
tific studies that describe the nature and existence of 
oceanic frents. Cromwell and Reid (1956) identified a frent 
by an abrupt kcrizontal density change at the sea surface. 
Roden (1976) Stated that a front should be considered to 
exist wherever any oceanic state variable (density, +tempera- 
ture cr salinity) reaches a relative maximum value in its 
Mori zenctal gradient. 

Piemterivartcn £Or studying fron=s and ocean prediction 
has keen stressed by Roden (1976), Elsberzry and Garwood 
(1979), and Niiler (1982). Elsberry and Garwood (1979) sre- 
cifically peint out antisubmarine warfare apvlications, 
fisheries management and a rescurce for ciimate research 
data, as viable and pertinent reasons for further study. 
Many important physical processes which are accompanied by 
attendant air-sea interactions occur in frontal tregions. 
Latent and sensible heat transfer, evaporative and preciri- 
tative processes, barcclinic currents, salt fluxes, mass and 
momentum transport, inertial-internal wave activity and mix- 


Pama l@ecanecceur. Oceanic frontal areas are areas of sound 
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speed changes and often are areas of increased biological 
Berl Vity. 


This thesis examines the time-depandent respons Od 


(DD 


huméerically- simulated oceanic front to local atmospheric 
forcirg. The model used in this work incorporates the forc- 
ings of wind stress and surface heating with a 20-level 
primitive equation ocean circulation model with embedded 


mixed layer dynamics. 


Be CESseRVALLICNAL STUDIES 

Cromwell and Reid (1956) described a simple case of an 
upper layer of warm, less dense water overlying a layer of 
Seow, Gencer water in the ccean. The front in this case is 
a density front formed by the horizontal temperature gradi- 
ent. Observaticns shewed surface convergence along the 
front, sinking motion in the immediate vicinity of the frent 
and divergence at lower depths. They observed that the wind 
affected the front by mixing the upper layer, reducing the 
horizcntai surface temperature gradient, and enhancing the 
¢*emperature jump at the base of the mixed layer. 

Voorhis and Hersey (1964) observed a wintertime thermal 
front in the Sargasso Sea. They found an average tempera- 


ture gradient of 1 deg cy/10 km, and a maximum in cne 


ie 





ioseance Cf 1.5 deg C/5 ka. This was accomparied by a 60 
cm/s jet flowing to the east alcng the edge of «he front. 
The depth of tke front corresponded to the mix2d layer 
ge pe hi The water on either side of the front wes found to 
b2 well-mixed. 

Katz (1969) continued the work of Voorhis and Hersey and 
derived teéemperature-salinity Geetelacsons for a ffcnt 
Observed in the Sargasso Sea in May. He found the frort to 
be a separaticn cf two distinct watar masses, with a tran- 
Siticn zcne cf only a few meters in the vertical. He 
deduced that water transpert along the front maintained the 
Shalp transiticn. Although not cbserved, Katz presumed mix- 
Pim MustecCCUL a> the interface, but the along-front current 
prevented any Euild-up of well-mixed waters in that area. 
He alse hypothesized that a vertical component to the cur- 
Bent MUSt Exist at the interface. Thus, the interface weuld 
be preserved by both horizontal and vertical motions, and 
the vertical velocity would account for the ebserved accuau- 
laticn of well-mixed waters below the front. 

Fang (1973) studied a front at the southern end of the 
Benguela Current. This was an intense front with a horizcn- 


tal temperature gradient cf 1.6 deg C/km and a temperature 
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emange 8OL 8-10 deg C across it. There was also local 
upwelling cccurring due tc an offshore transport. AS might 
be expected from Ekman theory, beng LOUNGES the -rone co 
intensify under southerly winds (Fkman transport being to 
the left of the wind in the Southern Hemisphere). Under a 
northwest cale, the front rapidly weakened. It rteintensi- 
fied es the wind steadied from the south. 

Regen has made significant and extensive contributicns 
more tre Literature on oceanic fronts, especially on «he 
large-scale fronts in the Pacific Ocean. In a series of 
papers, he has cbhbserved and eéxamined she Pacific subarctic 
meence (Roden 12975, 1977), the subtropical front {Roden 1974, 
1975, 1980), che doldrum front (Roden 1974, 1975), and «ke 
SUMakctac~subtrepical transition zone (Roden 1971, =1972). 
feecmce97 1 barez, Reden concentrated on the subarctic-sub- 
tropical transition zone and found that the wind stress 
Bavecean antcrtan= cele in determining the Location of heat 
and sait flux divergence zones. Roden (1972) focused on the 
temperature and salinity fronts located at the boundaries cf 
mye tlansa<Lcr zcne; the subarctic front to the north and 
the subtropical front (which is also a density tront) to the 


south. maemGrilGiis ct beth sronts were found to be related 
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to the wind stress, geostrophic flow fields and heat and 
salt flux divergence. Tke doldrum front was included in a 
pelayeGQt the subtropical front (Roden 1974). Further work 
en all three frents was cenducted by Roden (1975) with ref- 
erence tc the wind and to the energy flux fields. He 
observed that the depth of the mixed layer changed dramati- 
eaeeyecezO=ss the subarctic front in winter and spring; ‘fron 
100 mon the north side tc 300 mon the south side. The 
Subtrepical frent is highly dependent on the wind and energy 
flux fields at the sea surface, as the front occurs in an 
@rtea Gr EXMAN “LaNnSport ccnvergence. The doldrum front is a 
Beomicmaue te a Salinity gradient only, with a baroclinic 
2astward-flcwing current at the surface and a faster west- 
watd-flowing undercurrent. 

Hermerne Sikaret::c front, Roden (1977) found atnospheric 
MInMecingmat hat 15, Wind Stress, radiative heat flux, precipi- 
taticn and turkulent energy fluxes) to be the primary cause 
of frontal mcvement, intensification and decay. Roden has 
shown repeatedly that the roles of wind stress, Exman <*rans- 
port and keating are inextricable from the dynamics of fron- 
togenesis. noden (1980) again examined the subtropical 


Geentel zene {within which are located a number of 


i> 





individuai fronts) and showed that Ekman transport led tc a 
concentraticn of the temperature and salinity gradients in 
the upper cc¢an. Thus, frontogenesis was due to horizcnta 
convergence ard confluence of the flow field. 

While Roden's work has dealt primarily with frontogene- 
Sis ard in ¢€xplaining the dynamics of severai identified 
MeeenerPacatic frents, his findings are pertinent to this 
thesis, 1.@€., local atmospheric forcing plays a major rele 


meeene dGymawics cf oceanic fronts. 


.. MCDEL STOCLIES 

Gérvine has concentrated on a modelling. approach to 
fronts. In cne of his earlier papers (1974), ne investi- 
gated the dynamics of small-scale fronts, ~~ OUuGheew =< One 
m-gadrG to atmosSpheric forcing. SUGh es Ont SsDaopagatenanccr 
the ambient water via the unbalanced horizontal ovressur2 
gradient. Because of minimal cross-frontal mixing, he found 
Beromgescenvergence and Sinking at the front. He later 
expanded this rodel to include the effects of wind stress 
(Garvine 1979a, 1979b). However, because cf the need to 
Maintain steady state in the nodel, he prescribed wind 
strees to te cof secondary importance relative to the hori- 
zontal pressure gradient. Buoyancy and energy budgets were 


added later te the model (Garvine 1980). 
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Kao (1980) investigated the Gulf Stream as a large-scal 


(D 


density frcnt in the upper oceen. The heart of his work was 
performing a scale analysis, which ravealed three length 
scales; an inertial or deformation length scale, a buoyancy 
length scale, anda diffusive length scale. Kao's conclu- 
Sion was that the front was maintained by the cross-Gulf 
Stream circulaticn, but his model included no wind stress or 
heating. 

Cushman-Rcisin's (1981) model was based on Roden's work 
on large-scale frontogenesis, and it does have some similar- 
ities to the model used in this thesis because it includes 
mixed layer dynamics. An important difference is in the 
scales. Cushman-Roisin investigates large-scale and long- 
SeEM Spin-up cif a frent in respense to a wind stress curi. 
The model inthis thesis is used to investigare local 
respense of a pre-existing front tc local winds and heating. 
Surface wind stress is a primary means of forcing in bcth 
models, but Cushman-Roisin applies a strong wind stress curl 
MaemrGriect a wind field that changes direction at the 
latitude cf the zonally-oriented front. We, on the other 


hand, spin up a horizontally uniform wind stress over six 


ts 


hours, after which it maintains a constant value. The 
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results obtained by Cushmaa-Roisin depend much upon the curl 
of the wind stress. Cushman-Roisin finds the mixed layer 
depth is at a minimum inthe middie of the front, while 
Ekman downwelling is at a maxinmun. This is probably due to 
a lack of any wind stress at the front itself, aithough the 
Sel 25 Large. 

Camerlenge (1982) studied the large-scale response of 
Upper ocean frents, spécifically the Pacific subarctic 
front. His model has a variable-spaced horizontal grid with 
highest rescluticn in the frontal area. Vartical mixing was 
not included in the model. ASsuming a pre-existing front, 
Camerlenge arrlied different wind stresses, and found that 
except in the case of the passage of a strong cycione, the 
effects of wind forcing were limited to the upper 150 ao. 
One case cf a uniform wind stress did show Ekman transport 


and a horizontal displacement of «he front. 
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A. GENERAL 

The response of a pre-existing oceanic front to lccal 
atmospheric fcrcing (wind stress and surface heating) over a 
time period of two days is examined. An embedded mixed lay- 
er-ocean circulation model is used, and two different fren- 
tal structures are treated. fhe first numerically-inserted 
front, hereafter referred to as Frent 1, is somewhat of an 
artifact in that the isotherms beccme horizontal at a rela- 
tively shaliow depth (see Figs. 1, 2, and 3 for the initial 
conditions), so that Front 1 is manifested only within the 
upper 60 fm. The temperature field is vertically stratified 
belew appreximately 70 mo. The front is about 19 km wide. 
The paximum hcrizontal temperature gradient at the surface 
foils 9 dsq C/1.6 km. The mixed layer depth is a numerical- 
ly-inserted appreximaticn to what the actual mixed layer 
depth would be. MECCOSstrern- ec dlonmg-t=Onte Velocl Ey lsSmcal= 
culated frem the initial temperature field. It has @ maxi- 
Mum speed cf 32 cm/s with its axis located at the surface on 


the warm side cf the frent. The across-front geostrophic 
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Melocity is zero. PLON casi S saeViaple representation of ea 
smaller mescscale upper layer front, and the modsi results 
Mole §6©De “important in this context. The second type of 
front, hereafter referred to as Front 2, is a Simulation of 


a more typical front, similar tc those reported by Reden 


(19/80) . The horizontal temperature gradient at the surface 
Zemeles deg C/3.2 kn. The front is not only evident at the 
surface, but also its effect on the temperature field 


extends throughcut the thermocline of the model (ses Figs. 
mao oand 6 for the initial conditions). The mixed layer 
depth again is a numerically-inserted approximation based on 
the shape cf the temperature field. I+ is felt taat adjust- 
Hence wiieh cccur shortly after the start of integration of 
Eaemepodel will compensate for any errors in the initial 
epeld. Ree cctra nic along-Eront velocity 1s caiculated 
from the initial temperature field. Its maxinum value is 
100 cm/s and the axis lies at the surface on the warm side 
eee the frent. The acrcess-front gecstrophic velocity is 
Ze EO. Woene both Erehts oriented in a right-hand coordinats 
system, z increases upwards, y increases to the right, and x 


increases out cf the page. 
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Eve BYEOTHESES 

iemocmexecctea that the direction of the wind relative 
BOecne OfLenecaticn of the front plays an i#mportant™rOle in 
Bmemrecacrsetcntal circulaticn and density structure. Ekman 
transport acrcess the front should cause upwelling and down- 
welling features near the surface manifestation of «he 
me Ont. Jchannessen, ¢t al (1977) observed that «he thermo- 
cline threugh the Maltese front shallowed in certain areas, 
which indicated upwelling, and also that the thermocline 
Spatially deepened after isotherms surfaced. In Johannes- 
sen's case, this was attributed to the presence of an eddy. 
It the water masses had been of Similar types, it would have 
been indicative cf downwelling. 

Weetmeeecgeard t¢ One-dimensional mixing, WNiilerts (1975) 
mixed layer mcdel used an impulsive wind stress and no heat- 
ing, and showed that initially the mixed layer deepened rtap- 
ad Vs Strong inertial mctions can be generated which car 
lead to a larce increase in available turbulent energy and 
thus to a rapid(within the first half pendulum-day) deepen- 
Bogect che mixed layer. Mellor and Durbin (1975) also used 
an impulsive wind stress of 2 cm ys (approximately 2 


dynes/cm?) and ne heating. They showed that the mixed layer 
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Figure 2. Frent 1 Initial Mixed Layer Depth. 
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ae 


Geepened and ‘that “he temperature jump also increased 
Siaghcly. In a two-dimensional model, De Szoecke (1980) 
examined the effects of wind stress only in his version of 
the Niiler (1975) model and found that in the case of 
constant wind stress, tre mixed layer deepened rapidly in 


the first half pendulum-day (t=Piyt), and zhen slowed. An 


ta 


additional effect of the wind stress is the possible augmen- 


( 


taticn of a gquasi-gecstrcphic frontal shear (in che along- 
Peaont a, recticn) which could increase vertical mixing as 
prescribed by the initial value of the Gradient Richardson 


Number: 


de 
rv: 
—— < 


pa 


bP 


Ri = 


NO 


—_ 
ov 


Z 


oO” | “do 








A surface heating-ccoling cycle can lead to stratifica~ 
ticn under light to mcderate winds which will tend to modu- 
late any mixing tendencies brought about by wind-gensrated 
turbulence. WNiiler (1975) showed that surface heating halts 
the deepening cf the mixed layer. We expect the same result 
tn our model. Melicr and Durbin (1975) investigated a cass 


of acombined constant wind stress (2 dynes/cm-) and 3 
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constant heating/cooling value (H = +.01 cm-K/s). FOr teh 


(0 


heating case they found a shallowing and laveling off of th 


(wD 


depth of the mixed layer and an ever-increasing t2mperature 
jump value at the base of the mixed layer. The same heating 
condition with a lower wind stress (1 dyne/cn) resulted in 
a Shallower mixed layer depth and «the establishnent cof a 
double thermocline, a remnant of the initial mixed iayer and 
the fermaticn cf a new shallower layer due to the influence 
of heating. 

ieee seexpecc=d that inertial clzrents actoss the front 
may play an impcrtant role in the across-frontal transport 
of heat, mass, And moméentun. Phasing of inertial wind- 
Merven CULEEntS in conjunction with vertical mixing may 
mesult if laceral entrainment or “capture"™ of cee (warm) 
Mecer by the adjacent cooler water during herizontal oscil- 


Matmens acccss tre f£rent fliane. 


ce HOD 

To test the above hypctheses and to establish the basic 
dependence cf frontal response to local atmospheric forcing, 
four combinaticns of surface boundary conditions are inves- 


tigated. These are listed in Table I. 


ii) 
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TABLE I 


Atmospheric Fercing used in Model Runs 


Tx (dynes/cn_) Cy(dynes/cm ) Q(deq C-cm/s) 
Bae #120 0 0 
Case ITI -1.0 0 0 
Case [Il +1.0 0 -.004 
Case IV -1.0 0 -.004 


iMiemcne Cases above, Tis the along-front (x-direction) ccn- 


x 


ponent of wind stress, Ty PS. Cee LOcS-fiOnen(y=-direcei Ce) 
compcnent cf wind stress, and heat flux Q is applied to the 
océan uniformly at the surface (heat flux is positive 
upwards, sc the negative value indicates warming of the 
ocean from the air). All four cases are 2xamined fer both 
meaonce | and Frent 2. 

Cases I and II examine the effects of wind stress parai- 
lel tc the frent only, which will eventually generate a 
steady Ekman transport perpendicular to the front. The wind 
stress of 1 dyneycné corresponds <o a wind speed of 7-8 m/s. 
In all cases the wind stress is spun-up from 0 at t=0 to its 


Maximum aksclute value (1.0 dynes/cn*) by hOUWS 6.) Flom) Zou 


iD 


Pee oehourmmuc 2t remains ccnstant. Such a shett sSpin-up im 


(D 


is nex unreascnakle, as winds can change on a shorter tin 
scale, fcr example, when an atmospheric cold front passés. 


Cases III and IV impose a constant surface heating on each 
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Beene twC fFreceding cases. RSeendscaved oy the sanleer 
one-dimensional mixed layer mcdels, we expect that the hea-- 
ing value will act be overwhelmed by the wind stress, ani 
that the effects of both dynamic and thermodyramic forcing 
Will ke evident. 

Wiemin= ct: al conditions for Front 1 are shown in Figs. 1, 
2, and 3. Warmer, less dense water resides on the laft, and 
cooler, denser water resides on the right. The initial con- 
diticns for the velocity are determined by the associated 
Hylecal denszty field, which is balanced by a geostrophic 
peeomg=—chont CULDeEnt directed out of the figure (u positivs) 
msee Fag, 3). The mixed layer depth varies from a unifcra 
fermen im the tat field tc a shallewing through the front. 
The model is run for a 48 hour time period for gach of the 
specizied atmcspheric forcings shown in Table I. 

In Case I we expect tc see water from the cold (tigat) 
side transported to the warm (left) side, upwelling cn the 
cold side at the front and dcwnwelling due to convective and 
turbulent mixing cn the warm side. The mixed layer depth 
should deepen, and the front should diffuse. In Case II we 


expect the warm water to be transported on the surface *o 


miemeGeid side, but convective mixing to be absent due to the 
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stable conditicn of the less dense water overriding the den- 
ser water. Upwelling shculd occur on the warm side at thea 
front, and only minor downwelling should occur on the cold 
Side due to tke effect cf wind stirring alone. In both 
Cases III and IV the heating should have a shallowing influ- 
ence on the mixed layer depth tending to offset deepenin 
and cther effects caused Ey turbulent mixing. 


Pemonececal conditions for Front 2 are shown in Figs. 4, 


5, and 6. In this configuraticn, the warmer, less dense 
water resides ecnthe fright ard the cooler, denser water 
resides on the left. The accompanying along-front gecs- 


trephic current is directed into the page ({u negative) with 
its maximum value at the surface. The mixed tLayer depth 
Zollews the ccntcurs of the upper cthermociine, sShallowin 

Seao=ss the £frent. The mcdel is initiated and integrated in 
came eS 2t was fcr Front 1. 

In Case I (wind stress positive, directed out of the 
figure), we expect the Ekman transport to be toward the 
ect ts The results for this case should be analagous to the 
expected results for Case II applied to Front 1 as previ- 
ously noted. Case II should be compared to Case I for Front 
1 and Cases III and IV should ke compared to Cases IV and 


iererecapectively cf Front 1. 


52 





meer mercial rericd £¢cr our fodel is 


7. = exe (AG) alte) beats: 


27 
1 a 
euecne fatitude cf 36 deg N. Considering that in this mcdel 
wind stress dces not attain its steady maximum value until 6 
hours have elapsed, we should find a maximum in inertial 


moticn scmewha* later than half an inertial period after 


Ele Zoro, at rethaps t=15 hours or eo. 


D. Pesent@ektCn CF THE MOCEL 
The primary documéntation for «he model is in Adamec, ex 


OP Saga ik Cocr— 


al (1981). In that paper, the mcdél is cast 
dinates +o investigate the cceanic respons¢ to a Stationary 
radially-symmetric hurricane. In our case, the model is in 
right-hand Cartesian coordinates and it is two-dimersioneal 
oy and Z). It is believed that this model is particularly 
well-suited tc the study cr upper ocean fronts because it 


incorporates tke Garwood (1977) bulk mixed layer model into 


the Haney (1980) prinitive equation model. The ocean is 
assumed hydrestatic ard incompressible, where density is 4 
function of temperature aione. The Coriolis parameter is 


constant and there are no fluxes of mass, momentum or heat 
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MemeNeas tO any bctndary except at the surface, where atmos- 
pHerac £Orcing is applied. RicmEIOrizon- al Grid spacing 2s 
125 m and the vertical grid spacing is variable over 20 lay- 
ers from 6 m at the top layer to 64 m at the bottom Layer. 
The demain size in y and zis 100 km by 500 m, respectiveiy. 
ween tne fEcnt pre-existirg as an initial condition, the 
initial u-ccmponent of velocity is in geéostrophic balance, 
mache ieatial v-cCOmponent is lnitialiy set *o zero. 
Drawing upen Adamec, Ste eeCiooi,) nighliights. of the 
model are outlined below. Table II lists the pertinent 


model variables and ccnstants used in the governing 
equations. 


The ccverning equations in their two-dimensional form are as 





folicws: 
au ve wot 4 tv - Lt) (1) 
by wv tt ro - BR aS - Ea) 2) 
oa - 267) Gy 
a - — = 0 (4) 
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TABLE If 


Medel Constants and Variables 


o< coefficient of thermal expansion 
AT herizcental eddy ccerductivity 
_ coefficient for heat 
Ay Heaezcutal eddy @2fiusicn 
coefficient for momentum 
Rm vertical eddy conductivity 
 ,coefficient for heat 
K vertical eddy diffusion 
M coefficient for momentun 
g latitude 
ie Ccriclis parameter 


g gravity 
Zo BereLence density of Seawater 
To reference temperature 


ct 


mcdel timestep 

herizental grid array size 

Memeurcal Grid @€rray size 

basin depth 

rFasin width 

x (along-front) ccmponent of 
velccit 

y (across-front) component of 
velocity 

Tempera cure 

mixed layer depth 


prey < ahd Cle 
‘ 
o 


~~ 
i 


O0=- 2 -¢¢ (5) 


c 


The Ecundary ccnditions are: 


> . 
WO, fe = Anjo = 0 Cy = 0 , 


w=Q0@220, -D (8) 
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en Ge ea lee (6) 


20ba 27 6 ece 
22 5E+5 cn“/s 
2545 cm“ys 
0.5 eme/s 
ee cn¢/s 

8. 55 E=57=5 
981.0 cm/s¢ 
1.0276 g/cm? 
5 degC 

390 s 
80 

20 

500 m 
100 km 


100 km (7) 





Mek? = as , “u'w' = i » rviws =0@2=0 
O fo 


—w'??’ = -u'w' = -v'w' =0 @z=-D (9) 
An eguaticn fcr the depth, h, of the well-mixed layer is 
derived by integrating the continuity equation and applying 
@eemrs2aqtdq iid kcundary condition (Eq. 8): 
+ = 
i ae ee (19) 

The vertical turbulent fluxes are parameterized in two 
dirferent mtanners depending upon where in the water column 
marey acre located. Below the mixed layer, the vertical 
fluxes are parameterized ty eddy viscosity and conductivity 


Goecrtficients: 


ea gee ae t ou 
yiwf = “hss 
viw! = -K,S° G4) 
wit! = -K noe 


Above «he mixed layer, turbulent mixing makes such a method 
which assumes Kur an d Km constant unrealistic and inaccurate. 
Thus «he Garwecd (1977) turbulence closure model is invoked 


using bulk turbulent kinetic energy equations. 
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The entrainment buoyancy flux is given by: 
ees sis 
nKgwit')_, = ~ Wwiw') = &)/h (12) 


where <w'w'> and <E> are the vertical and total compconencs 
respectively cf the turbulent kinetic energy. The mixed 
layer total turbulent kinetic energy equation and the verti- 


cal component of turbulent kinetic energy Equation are: 


2 PS - nu? - agh(wT*)_,/(2Ri") = ( a)? + fh)G@) (13) 


2 = agh((W'T*)_, - (WT) )/2 
+ () = 3qutw')) G)* 
z ((3)* eyes) (14) 


aerenceute = ex 2 ieomenee fe iC-1OU Volocicy and 


Ra 


Ri*= 0 is the bulk Richardson Number. 
Aus +ave 


The downward fluxes of momentum -(u'w') and -(v'w') are con- 


° ES Ena terS ea: Au 
(15) 

V 

eee sel (wt) (BE) 
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where w, is tke entrainmert velocity: 


e 


ee, 


AT (16) 


The entrainment heat flux COR is the lower boundary 
condizion for the mixed layer temperature profile, and a new 
mixed layer depth and new momentum and temperature profiles 
are calculated. The mixed layer depth is independent cf the 
model level derths. If the mixed layer shallows rathar than 
deepens (as wculd occur if there were Rec teas Cae ee ee | 
turbulent energy to transport heat down to the depth of the 
existing mixed layer), then the depth where (w'T')— 0 
becomes the new mixed layer depth. During mixed lay¢r shal- 
lowing, heat, momentum and potential energy are conserved. 
To prevent the levels beneath the mixed layer from becoaing 
unstable in the numerical frocess, a dynamic stability con- 
diticn is impos2d such that the Gradient Richardson Number 
2S always greater than 1/4. 

The embedding of the mixed layer model within the gen- 
éral circulaticn model réguires communication between the 
two. The general circulation model is the dynamic portion 


which calculates changes of u, v, and Tat each depth level 
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due to advective and difftsive processes. These valuss ar 
then acted upen Ey the mixed layer model, which in turn cal- 
culates the changes due tc surface fluxes and entrainment 
Mixing, and calculates the new mixed layer depth. These 


values are then used Fy tke general circulation model. 
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III. MODZI RESULT 
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A. ANALYSIS AND DISCUSSICN 
Miemsoudiecne LEOnt 15 anitially in geostrophic balance, 


the impositicn cf atmospheric forcing causes the front to 


undergo a variety of adjustments. In Case If, 4 positive 
Wind stress cf 1 dyne/ene is yea pelied . The front and the 


M—cOmponent ficid undergo diffusicn, thereby giving rise te 
an ageostrophic u-component, which in turn creates a v-ccn- 


ponent. As diffusion occurs through mixing, and 2s the sur- 


face expressicn of the front is advected to the left in 
respcnse tc tke wind stress, the along-front jet decreases 
in magnitude as well as in horizontal extent. A small 


countercurrent is present at hour 24 dué to a reversal in 
the siope of the isctherms in the vicinity of the front. 
This effect disappears by hour 36, but then reappears at 
hour 48. 

The pase cf the front cn *he right side deepens (as evi- 
dent by the isotherms beccming vertical to greater depths) 
eeameapeuc 3S0)n inaztialiy to abcut 50 m as a result of the 


effect of turbulent mixing and the creation of ageostrophic 


4Q 








velocities while undergoing dynamic adjustment (se2 Fig. 7) 


} 


An ageostrephic u-component is evident in Pig. 8 in the area 
of y=53 xm to 61 km. The mixed layer depth increases as the 
upward slicfing thermccline erodes within the front because 
Mieeteirs transported across the front, right to left. The 
mixing of the cocler water with the warmer water causes the 
mixed layer depth to deepen on the left (warm) side. The 
axis of maximum v-valiues lies just above the bottom of the 
mixed layer near y=53 kno. If the v-profile is @xamined 
(Fig. 9), it is seen that dv/fdy > 0 on the tight (cold) side 
of the frontal axis, and dv/fdy < 0 on the left side. The 
gradients cf v are much stronger on the left side of the 
front than they ere on the right. The maximum v-velocity is 
omic |O6NeCM/S. FrOM continuity, dw/dz < 0 (downwelling) on 
meme E2ght side cf the axis and yee > 0 (upwelling) ina 


thin layer on the left side that closely follows the mixed 


layer depth. It is noted also that the mixed Layer decch 
has deepened at the right side of the front. [his Gone ven— 
ancn is oresent throughout the model integration. Though 


Such a conditicr was not anticipated prior to the running of 
the model (we expected an upwelling in the mixed layer on 


the right side as a result of the mass transport over the 
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front towards the left) 


believed tc be physically 


this increase in the mixed layer depth occurs 


non-staticnary adjustment 


also occurs in Case II w 


stress is applied. Becat 
+o the 
challcwer 


mixed layer tis 


field plocks 


water enccunters this "obstacle", 


Side and thus fcrces a deéfression in 


explanation must ke temper 


are instantanecus snafshot 


vertical integral of 


<ne £low fren 


this depression in deéeptn is 


g 


ce 
tow 


realistic. It is suspected that 


in response to 


tO the a@tne@eonaeric forcing, as it 


here the oppositely-directed wind 


se the mixed layer depth resronds 


pRanspeGe,. the Lact the 


- 
— 


Mw teens che ofcoen= chen in “the fag 


. 


penetrating to depth. As th 


( 


Pemmoacksaip son tie agit 


the mixed layer. This 


ed by the fact that these figures 


¢ in time, and can be deceiving in 


attempting te descrite processes that have been evolving 
over a length of time. Another possibility is that some 
type cr internal standing wave has been excited which is 
causing fluctuations in the mixed layer depth. 

The influence of the ageostrophic u-velocity, the core 


eo «~he Vrweloci~ 
depression uf the mixed la 
is felt 


welling are evident 


Pana cage Cf the front. 


lying along the weli-mixed 


at depths greater than 


at these depths 


layer and ths 
yer on the right side of the front 
100 m. Upweiling and downr- 


CRG. 7) beneath the 
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Figure 7. Frent 1 Case I 12-Hour Temperature. Contour 


interval is 0.z deg C. 
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Figure 9. Front 1 Case I 12-four V-Veiocity. Contour 
interval is 2 cm/s. Dashed lines indicate 
negative values. 
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Interesting patterns in the v profiles evolve a+ hours 


(n 


24 and 48. At 24 hours a small positive v exists. Teese 


Semeraty tO the expected Ekman trensport direction in th 


@ 


surface layer, and is protably the result of an oscillation 
created by the inertial mction. Furthermore, at both 24 and 
48 heurs, two cores of a positive v-component are present 
symmetrically centered about the base of the mixed layer on 
either side cf the front. These are not present in the 
front, but form outside the frontal area and extend =o the 
far field where the mixed layec is reiatively stable. 

The across-front velocity component 1s in conscnance 
With the directicrn of the net Ekman transport at all times 
excerpt hour 24. As noted, this velocity 1S at a maximum at 
Smeewbase Of the mixead layer rather than at the surface. i 
is nen-zerce Leneath the gixed layer in the middle of the 
Beont as well. Mpc ~w etc appedrs that interfacial mixing 
does occur in this case and it may be contributing to the 


vertical mixing cf heat in the thermocline at this location. 


Die Suman e—-—- CASE I 
Under the influence of positive wind stress and no heat- 
ing, the frent has diffrsed and has been advected in the 


direction of the surface Ekman transport. Because the frent 
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is unsteady under an arrlied wind stress, ageostrcphic 
velocities are created. Faas ©@f£ zhe seSponse is fricticnal 
due to the turbulent boundary layer or nixed Layer oro- 
cesses, and part is inertial due to the applied wind stress 
PouncaLy ccnditicn. A neticeable deepening of the mixed 
Mawes OGCUrS at the front eon the right (cold) side, and the 
horizontal variability in mixed layer depth and temperature 
are reduced significantly under the influence of both verti- 


Gal and excess-frontal mixing. 


Qi] 








A. ANALYSIS AND LCISCUSSICN 

In this cas¢, a negative wind stress (directed into the 
page) is applied. Wind stirring must overcome static sita- 
bility as the warmezx water will override «he cooler water. 
It will be seen that the wind stress cf 1 dyne/en? £5 SuUtfii-— 
Cclent to mix the warmer water which is transported toward 
the ccoler water. There is no evidence of stratification of 
shallow waters. Recall from Case I that, as the denser 
water from the right side was being transported +o the left, 
2t had to overceme static stability in order to upvwell and 
override the front, and that the mixed layer depth profiles 
posed a vertical barrier tc the transpert. In this cas, 
the base cf the mixed layer is at an incline in the direc- 
tion of transport and static stability need not pe overccme 
in the transfert process. This leads to the expectation 
that cross-frorntal mixing may be minimal in this casé. 

Diffusion in the temperature field again is evident over 
the 4&8 hours, but this tire it occurs primarily at the edges 


Greene fEeOont. A strong horizontal temperature gradient of 
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0.6 deg Cy1.2 km in the middle of the front is maintained up 
through 48 hcurs (see Fig. 10 as representative of this) 
which was not present in Case I. This lends credence to the 
hypothesis that cross-frontal mixing is minimal. 

Moerceable Mixing occurs on the cold side of the fecn<. 
The isotherms are vertical to greater depths than they were 
in Case I (see Fig. 10). The turbulent mixing alone not 
Smy CVGECCMSs static stability but is also able to mix «c a 
relatively greater depth cn the right side of the front. 
The water that is being transported across the front is 
apparently immediately mixed by wind stirring as no plume of 
the warm, less dense water is present on the right side. 

The alcng-front current again forms an ageostrovhic ccn- 
ponent as it urdéergoes adjustment, and is especially evident 
See ne siagnt Side of the front around y=55 kn (Fig. 11) 
where it exists to a greater extent both horizontally and 
vertically than was seen in Case I. The core of the alcng- 
front jet dces decrease in speed and diffuse in time. BY 
hour 36 the maximum value cf u has been reduced to slightly 
more than 8 cmys. The creation cf this agecstrophic veloc- 
ity accompanied Ly a v-compenent leads to strong downwelling 


around y=55 km underneath the front (Fig. 10) and remains 
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Figure 10. Front 1 Case II 36-Hour Temperacure. COn = Cus 
interval is 0.2 deg C. 





over the 48 hcur integration. A strong core o£ negative 
across-front velecity of cver -6 cm/s is created in «his 
ateéa (Fig. 12), and there is an attendant sharp upwelling 
in the isctherms at 75 m and wave-like fluctuations present 
to depths beycnd 175 m In the areas where dv/fdy > 0 (dw/dz 
< 0), downwelling is seen in the front, and where dv/dy < 0 
(dw/dz > 0), urwelling is seen. In the surface layer, the 
v-component is positive and in the direction of Ekman trans- 
port. It is ccncentrated in a jet in the area of the tight- 
est horizontal *emperature gradient over the shallow peak in 
the mixed layer depth. All of these tendencies and trends 
persist over the 48 hour integration. 

A deerening in the mixed layer occurs to the right of 
the front near y=55 km and progresses outward in time. The 
Shallew peak in the mixed layer depth within the front is 
mostly preserved, though somé erosion and mixing obviously 
has cccurred (ccmpare Figs. 13 and 14). PRES “Gee ner .sule- 
Potts tne hypothesis that minimal cross-frontal mixing 
occurs. Examination cf the corresponding v-fields aise show 


[ieee lesGcEces-trontal tixing. 
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Bie Soe s -- CASE IT 

In this case of a forcing of a negative wind stress and 
no heating, the temperature structure of the front retains a 
strong horizental gradient over the 48-hour integration. 
Pecmaicr on Cf the front occurs, but it is concentrated at the 
boundaries of the front. kind stirring overcomes the static 


transported water once it crosses the fren. 


rh 


a2 


j-- 
ct 


stab vc 
The creaticrn under adjustment of a large ageostrophic u-va- 
Moczty amd an accompanying v-velocity give rise to strong 


upwelling and dcwawelling features within «h2 front. The 


shallew peak cf the mixed layer depth within the front is 


j4- 


iD 
3 
I 
(0 


preserved cver tin fact that this peak is maintained, 


(+ 


moe herizcntal emperature gradient retains a tight struc- 
mr=e, and the vy~-field 1s mainly constrained to the agecs- 


trophic areas and surface layer, suggests that interfacial 


Semecrcss-trcntal mixing 1s Mininal. 
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Pigure 12. Front 1 Case II 12-Hour V-Velocity. Contour 
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Figure 13. Frent 1 Case II 24-Hour Mixed Layer Depth. 
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V. MCDEL RESULIS FOR FRONT 1 CASE III 


A. WHALYSIS AND CISCUSSICN 

A unifcrm surface heating is applied with the identical 
Wind stress cf Case f. The magnitude of each forcing is 
such that influences frem both are expected to be mani- 
fested. In examining Figs. 15 and 16, the obvious effect of 
heating is tke adjustment of the temperature structure on 
Seenece Sldae Cr thre» front. At first examination, it appears 
that the heating is having the opposite effect of what is 
expected — <that is, it appears to cool the mixed layer 
rather than warm it. In reality, the heat is carried dewn 
to depths cf 50 mand warms the water there. Cooling occurs 
between 32 and 5C m depth, which is readily seen in Fig. 17. 
Mixing may create a dynamic instability in the thermocline, 
and vertical diffusion of heat by Km #42 T /dy ¢ may cause the 
intermediate water to be ccoled and the deep mixed layer to 
be warmed. This basic structure remains over the 48-hcur 
Mee GLatacnh. The frent undergoes diffusion and is advected 
to the left, as cccurred in Case I. The front is manifested 
tO a lesser deépth in acccerdance with the heating process - 
downwéelling and mixing have been moderated by the heating. 


Though some wave-like mction is seen in the temperature 


oF 





peetd beneath the mixed layer, the magnitude is not ae great 
as in the case of wind stress alone. The mixed layer d¢oth 
has shallcwed tc about 43 m under heating alone in the areas 
outside the frent (see Fig. 16). This shallowing occurs 
because the atmrcepheric fercing (heating and wind) is input 
as a body fcrce cver the entire mixed layer. The heating 
has the effect of making *he mixed layer depth more uniforn. 
The marked deefering in tke mixed layer depth that occurred 
in Cases i and II at the right (cold) side of the front is 
not present in this case. The deepening of the mixed layer 
mertectS =he u-field. The agecstrophic u-component does not 
develop to as qreéeat an extent as it did in Casa I. The gen- 
eral along-frent current eppears much as it did in Case I, 
Pe eeedate t= (vertical extent is also conserained by the 
Shallower mixed layer depth. In both casés, maximum speeds 
agree cover the 48-hour integration, and they r2 advected 
with the front as it respcends to the mass transport. 

The v-field follows the same trends as i+ did under wind 
Stress alcne. At 12 hotrs, vis directed to the left as 
expected following the transport, end at 24 hours the v-ve- 
locity reverses and establishes the layers of concentrated 


positive v-values centered symmetrically abcut the mixed 
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Moerion= | Case LII - f vs z Profil 
taken at y=35 km for t=0 and ~=1 





Meyereaepth. At 26 khcurs, v returns to a negative ve 


}— 


ué 42g 
gep4O HOUrS, though still negative, vis much decreased in 


magnitude. Again, cells cf positive v-velocities ar 


(iD 


pres- 
ent at about the mixed layer depth. The across-front velcc- 
ity does fenetrate the front and the shallow peak of the 
mixed layer (s¢e Fig. 18) which is indicative of cross-fren- 


tal mlxing. The shallow séction of mixed layer itself is 


ct 


eroded and smccthes cut cver the 48-hour integratior as i 


Geceetnecase I (Fig. 19). 


B. cur, -- CASE IIT 
The effects cf bcth wird stress and heating are cleat in 
this case. Tke wind stress induces a mass transport to the 


right, moving ccoler, denser water into warmer, less dense 


Wa =CL. This water not cnly uapwells and moves over the 
front, but also moves acrcss the frontal interface. At 24 
and 48 hours, the v-field appears <*o be influenced by 


inertial oscillations with a reversal in directicn occurring 
in the surface layer at 24 hours and a hnear-reversal at 48 
NOUES. Cores cf positive velocities are established cen- 
tered about the mixed layer depth outside of tne front. The 


front bozh diffuses and mcves to the left in time. 
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H at hour 48 
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Heating shallews the mixed layer depth as expected, and 
thereby restricts the vertical extent of the velocix 
E2elds. Heating also prevents the immediate deevening of 
the mixed layézxr on the right side of the front, which was 
also seen in Eeth Cases I and II. The heat manifests itself 
at the initial mixed layer depth where the isotherms become 


more stratified. 
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VI. MODEL RESULTS FOR FRONT 1 CASE IV 


A. Ans LS LS AND PESCUsS1e N 

Case IV ccmbines a negative wind stress and e2 uniforn 
surface heating. The mixed layer shallows noticeably and 
heat is carried +o the bcttom of the mixed layer 2t depths 
Seeoo co J0em aS in Case III {see Figs. 20 and 21). The 
heating dces have an interesting effect when compared with 
the wind-stress-cnly forcing of Case MII. Pe ashiae wy da ous 
the isotherms in the middie of the front is not evident here 
as it was in Case II. H¢rLZOheel Gauls won of she tenoera— 
ture structure cccurs more rapidly and attains @ gradient of 
1.6 deg C75.6 km by hour 12 (see Fig. 20). Give 2 ong -ssenc 


velocity is ligited in its vertical development on the right 


iD 
rh 


saae Or th Bete Veal Y=S? knee. 22) and follows biech 
Sve iGerzZonetal ard vertical diffusion tendencies of the’ den- 


Saeiy LEGnet. The across-front velocity has a narrow, verti- 


cally-oriented cera of negative values with a maxinum of -6 


th 


cm/S cn the right side cf the ront between 25 and 50 o 
d= pemeasean Cass IT. Since this maximum of v-velocity 


exists in an area where there is little u-velocity, its 
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existence cannot be attributed to a response to the agecs- 
trophic u-ccmecnent. The positicn of this v-velocity corre- 
lates well with the upwelling/downwelling pattern evident in 
the isotherms in Fig. 20, and also with wave patterns in tha 
isotherms over the entire depth extent in the right half of 
ds 236 We suspect that such a large area of v-component 
created below the mixed layer is a response tc the warming 
of the water at the initial mixed layer depth which has 


i 


9) 


shoaled, re-oriented the isotherms and created vertic 
velocity. Furthermore, wind-generated mixing forces verti- 
cal motions in the water mass tc a great depth, inducing a 
v-corfronent there. 


i'peerec case, the front diftuses horizontally more rap- 


th 


idly «han under wind stress forcing alone. The = Che 
spreads like a warm, iess-derse plume overriding the cool, 


denser water (ccmpare Fig. 24 t9 Fig. 20). Ehie 2 long-f rene 


{Aa 


GWEGSEoeSHAastS With the front. The mixed layer loses its 
Seoowesect2 On Within the front and smoothes out, which did 
not happen in Case II. The peak in the mixed layer was 
maintained much lcnger. With heating, +here is less of a 


temperature difference and less of a shallowing from the 


far-fisid mixed layers to the shallow mid-frontal mixed 
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Figure 20. Front 1. Case IV 12-Hour Temperature. Contour 


interval is 0.Z deg C. 
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Figure 21. Frent 1 Case IV 12-Hour Mixed Layer Depth. 
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Figure 22. fFrent 1_ Case IV 12-Hour U-Veslocity. Contour 
interval is 4 cm/s. 
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V at hour 12 
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Haver, which frebkably acccunts for the quicker breakdcwn of 
the shallow section of the mixed layer. 

The v-field in the surface layer is advected to the 
right, with a maximum value located at the surface within 
mene front. By heur 24, strong horizontal cells cf a v-ccn- 
ponent have fcrmed about the mixed layer depth. In this 
case, the values are negative, whereas in Cases I and III 
they were positive. (Case II never formed them.). These 
cells disappear at hour 36 and reappear at hour 43, as they 
did in Cases randee bi. The mechanism that explains the 
appearance cf negative v values is unclear. 

The trends to 48 hours show further diffusion and advec- 
tion of the warm frontal plume to the fright. There is 
SOGtesnemece the herizental G@aistribution of mixed layer 


depth. The surface front patterns of u- and v-velocity are 


all advected te the right. 
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T at hour 36 
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Peguteecweeerren: 1 Case IV 36-Hcur Temperature. Contour 
Pimereval is 0.2 deq C. 
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B. Surah == CASE IV 

Applicaticn cf a negative wind stress and surface heat- 
ing results in the front moving @s a warm cverriding plume 
to the right. The mixed layer initially shsllows due to the 
heating, and cver the 48 hour integraticn breaks down and 
SMOcCthes CULT in the herizecntal. TRE along-front velocity is 
advected and diffused in the mixed layer. frien pei Cranes 
memoce- rae: vVelccrty shifts to the right in the surface 
layer, with its maximum value centered over the shallowest 
portion of the mixed layer. A v-ccmponent develops to quite 
Beiecacmadceeh Ci the Lignt Side Of Erecnt due to motion gén- 


erated by heating and turktulent mixing. Celis of neqativ 


Uy) 


v-comrponent are fermed at hours 24 and 48 which are centered 


about the mixed layer depth. 
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VII. MODEL RESULTS FOR FRONT 2 CASE I 





A. ANALYSIS AND CISCUSSICN 

Dnementadt CORdi tions for Front 2 are shown in Figs. 4, 
wee and 6. The front extends throughout the depth of the 
basin and is criented such that warmer, less dense water is 
on the richt and cocler, denser water is on the left. Many 
of the dynamical features of Front 1 are expected to ke 
repeated fcr Frent 2. ENBCOMPa Gingmh cont 2 ©O Fromt 1, the 
Pome zcncal axes are oppositely-directed <¢o0 the frontai ori- 
entation. The thermocline slopes upward in the positive 
Mecmtectacn fOr Front 1, tut downward for Front 2. 

The wind stress cf Case I applied to Frent 2 creates a 
Pane Pere Girected to the lert. Warmer, less dense water 
moves towards the cooler, denser water. iMenet = S Gass, dl t= 
fusicn in the temperature field is present for the duration 
of the integration. The horizontal temperatures gradient 
remains large (1.0 deg C/1.2 km) in the middie of the front 
(aac neme Majesty OL the diffusion occurring at the right 
boundary. Meng On the cool (left) side of “xhe £rcnt 


extends £rem the surface to 135 m-which is below the initial 


aS 





depth of the mixed layer. Resuits of downwelling are seen 
at the left frontal boundary, near y=42 kn. Apparently, 
there is no chance fer a stable, overriding warm plume to 
develop before it is mixed into the ambient cool water pecl. 
These phercmena in the temperature field were also seen in 
Sase [2 fer Frent 1. Much more dramatic oscillations occur 
in the iscthéerms, primarily on the right side between 53 kao 
and 70 kn, where the entire field is inclined. They area 
also visible cn the left side superimposed upon an inclina- 
ticn of <he entire field sloping downward to the right (Pig. 
Bh. The isctherms show a definite trend towards becoming 
more uniternly irclined beneath the mixed layer across the 
entire horizcntal expanse, with the steepiy-inclined field 
@eeeene Eecht Leccming more horizgortal. Downwelling causes 
the lef+ side tec slore dcwnward to the right and upwelling 
lifts the isctherms on the right side. 

The mixed layer depth undergces major changes over the 
entire domain. Cn the left side, it deepens to where the 
depth in the far field is approximately 125 m (Fig. 26). 
Since this deepening occurs by 3 hours of integration, such 
a drastic deepening probaktly resulted from the initial mixed 


layer depth being inserted at too shallow a depth in this 
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area. As in Frent 1 Case II, a downward bulge in «he nixed 
iayereextenasng tO 155 m Cecirs at the left frontal boundary 
at y=43 km (Fig. 26). In this case, the downward bulge does 


not ccrresrond to the side on which the waters would «end to 


Bede k US" Pxrlcr to Crossing over the front as in Front 1 
Case i. As strong mixing is implied by the vertical orien- 
taticn of the isotherms at that location, the bulge is 


believed to ke associated with a cooling of the waters in 
that area as a result of enhanced vertical mixing caused by 
the ageostrophic along-front velccity shear (Fida. 27) % 
Within the frent, a shallcw peak in the mixed layer depth is 
Mainteined at the left bcundary, though over time it does 
deepen sligh«=ly and erode through mixing. There is much 
variability in the mixed layer depths on the right side cf 
enc LEOnt, but it appears *hat the general trend is for 
Geepcpingecs cccur in the far field beyond y=55 km (Fig. 
26). 

Along with the isotherms, MMe malomg—LLOnt  U-Vve Locicy 
diffuses ard décreases its maximum speed over time. As 
already seen, an ageostrophic component forms on the lef- 
frontal bcundary and stimulates mixing there. The remnants 


of a tight herizcental velccity gradient are seen after 48 


78 





Depth (m) 


H at hour 36 
0.0 


25.0 
50.0 
75.0 
100.0 
125.0 ™ 
150.0 
175.0 
200.0 


225.0- . —— 
25.0 35.0 45.0 55.0 65.0 75.0 


Y-distance (km) 


Figure 26. Frent 2 Case I 36-Hour Mixed Layer Depth. 
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Figure 27. 


80 





foleeyeonewgn ayecle cr ne advecticn of the along-front cur- 


rent occurs in response to the wind-generated +*ransport. 
Wave-like cscillations aprear to propagate along the deeper 
extremities cf the u-fielid, but no phase relation between 
isotachs cr with cther fields is readily apparent. 

The acress-front v-velocity does correlate with the 
Mire cr.on Cl transport ur threugh hour 9. An oscillatory 
Barmesvetr @xac Ss ir the w fieid over the 48-hour integration. 
Ac hour 12, the entire v-field is positive, but this may be 
a response to the agecstrophic u-velocity, and it is streng 
enough to cvercceme the acrcess-frent transport velocity. At 
PeePours, VY is in the direction c£ transport and a maximun 
ot -20 cm/s appeers at the surface within the front over the 
Shallcwest secticn of mixed layer depths at y=43 km (Fig. 
Zo}, and temains through hour 36. This is reassuring as it 
was expected that the inherent static stability of the water 
masses would minimize deep cross-frontal mixing, similar to 
what was seen in Front 1 Casée If. The isotachs very defi- 
nicely terminate at the left boundary of the front, but they 
also exist cbeneath the mixed laver and co cross the front 
beneath the mixed layer. At hour 48, v is positive again, 


perhaps due tc inertial mctions. teas proposed (Sheugh it 
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fomececrtainiy fot as clear as it was for Front 1) eet 
cross-frontal mixing, PeelGdwiniy pieced, OCCUrS here <tc a 
Saeates sextent then st did for Front 1. The majority of the 
Mass transport appears to exist in the mixed layer and 
passes over, rather than through, the shallow mid-frontal 
depth of the mixed layer. The reasons for the variaticns 
between the fronts may be due to the greater u-componéent 


(both in speed and in areal coverage) present with Front 2. 


B. SUumaey == CASE I 

Several ccnuplexities have been incorporated into the 
modeling of Front 2. First, the isotherms slope dewnward to 
ee tLegncs cver the entire right half of the basin. Hence, 
SemimonrenO: Chiy exists at the surface but is "connected" 
through the depth of the Easin. This establishes an initial 
u-velccity which exists over the same area. Adjustment to 
the wind stress may not be confined solely to the upper few 
meters, as dynaric respenses may be transmitted to depths 
beyond those directly affected in the surface layer. Also, 
the mixed layer depth is artificially inserted és an initial 
Somaettco,eand Tay net be the mest perfect fit to the temp- 
erature field. Though this mixed layer configuration may 
exaggerate cr suppress actual features, it is nevertheless 


felt to be a reasonable approximation of an actuai front. 
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The responses are much more complicated in structure in 
Front 2 than they were in Front 1. Bayond 12 hours, the 


front shows little advective response to the mass transport. 
One reason this cccurs is Fecause the front exists over such 
a great depth. The criginal thermal structure and temnpera- 
ture gradient are maintained in the middle of “he frort, 


ae bHOrLizental dzifusicn cecurring primarily at the right 


boundary. Mixing in the surface layer occurs as the isot- 
herms beccme vertical to qreater depths over time. At thea 


fi 


left koundary cf the front, downwelling is quite evident an 
the waters mix easily due tc the creation of an ageostrophic 
u-component which enhances the wind stirring. Wave-like 
osciliaticns are evident in the temperature field beneacéh 
the mixed layer, and the entire deep temperature fi¢id 
becomes mere uniform in slcpe over the 48-hour integration. 
The mixed layer depth must be considered in three sepa- 
Patemcegimes = left side, frontal, and right side. On the 
left, the mixed layer depth increases by almest 50 m at the 


frontal beundary, where mixing is enhanced due to the influ- 


ence of the acecstrophic current. The mixed layer bulges 
downward, then rises sharply to its shallow peak at thse 
Ape(op dle This shallow section remains throughout the 48 hour 
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integration though it dces deepen slightly due to mixing. 


On the right side, the mixed layer deepens and becomes mor 


TP 


horizontal thcugh there are strong fluctuations over inter- 
vals of 4 km and iess. 

Iwo primary trends aré observed in the along-front u-ve- 
locitv: 1) tke u-component diffuses in time, still retain- 
ing its basic shape but decreasing in speed and horizontal 
gradient; 2) an ageostrophic ccmponent is evident, prima- 
rily extending cutwards frem the left boundary of the front 
and which is tkelieved to "drive" the downward bulge in the 
mixed layer there. 

THevacross-front welccity ccmponent 1S ccncentrated at 
the surface in the direction of transport between 0 and 9 
hours and at 24 and 36 hours. The front poses no barrier to 
the v-velccity at 12, 24, or 48 hontrs. The 36-hour fields 
apeue=mpe Crly cres which have essentially no cross-frent 
Demiceavymcen vOrer. Other than im che surface layer. Because 


the initial herizontal temperature gradient is preserved in 


© 
rh 


the middle cf the front, and because the shallow portion 
the mixed layer depth is still evident after 48 hours, = sug- 


ese tnawele st le Cross-frcntal mixing occurs. 
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VIII. MODEL RESULTS FOR FRONT 2 CASE II 


A. ANALYSIS AND CISCUSSICN 


With the wind stress cppesite from that of Case I, «the 


i | 


CoIcerntMeaeetlanspcre is left to right, and cooler, denss 
water will be transported «oward warmer, less denss water. 
There are Similarities here tc the results of Front 1 Case 
is The front becomes much more diffuse at the surface and 
the horizental temperature gradient is not preserved as it 
was in the preceding case (see Fig. 29). It has decreased 
Eomiess ageg C/6.4 km. The hecrizontal field on the left side 
reorients Ey inclining down to the right, and the inclined 
field on the right becomes more horizontal. Trae. total 
effect is tc make the terperature field more level beneath 


the mixed layer. Wave-liike fluctuations are also seen in 


cr 


h 


(D 


isotherms (Fig. cove Vertical mixing can be readily 
seen «o have cccurred within the front by the steep slope of 
the isotherms. 

The maxed layer depth immediately adjusts in the letx 
preooocmce ereunrd 125 mas it did in Case I, but a much 


greater dcwnward bulge in the mixed layer forms at the left 
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Figure 29. Frent 2 Case II 36-Hour Tempstacure. Contour 
interval is 0.5 deg C. 





boundary cf tke front around y=42 km and remains over the 48 
Hows sntegration (Fig. 30). (Mietoesesult ats Consisten* with 
the interpretaticn preposed for Front 1 Case I, where tne 
Shallower mixed layer within the front acted as 4 barrier to 


mass transport, causing the transported water to "back up". 


tt 


The mixed layer peak within the front deepens and erodes. 
The right side far field is again variable but shows the 
same tendencies to deepen. 

The agecstrophic u-velccity (not shown) is vary similar 
to the iscetach pattern which evolved for Case I. AVieecre Or 
maximum velocity is advected to the right which correlates 
Wesieotemacvec-o2cn of the front. The biggest differance is 


that the isotachs appear at a greater depth just inside the 


Pettey L=encval Ecundary. 


fs 
W 


Miemexgected directicn of the across-front velocity 
PE eecOmetenu, ={he directicn of the nezt Ekman transport. As 
someaseuwes, a wave oscillation is evident. Atei2 andra 
moms =n LS Clearly so, with the velocity field existent 
along the entire section ci the front from the base of the 
PIOmceateawetnd 150 m ip to the surface (Fig. 31). The max- 
imum velocity cccurs ata depth bernsath the shaliow mid- 


meonedemeseGtscn Cf the mixed layer, which s¢emns <0 
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Figure 30. Frent 2 Case II 48-Hour Mixed Layer Depth. 
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Sole sesieerteo Cross-frontal mixing, erosion of the nixed 
Maygec aGepth watchin the frent, and weakening of the frental 


horizontal temperature gradient. At 24 hours, the v-velcc- 


j de 
{-4- 


+y field is reversed, and negative velocities predominate 
in the frert and beneath the mixed layer. Inertial motions 


probably account for this reversal. 


B. Ssdmmant ~—- CASE II 

The temperature structure of the front shows vertical 
mixing and hcerizontal diffusion cccurring over the 48-hour 
integration. The horizontal temperature gradient is 
teomedscqme aque to =he effect of convective mixing cf the 
cooler water which 1S ‘transported to the warmer water by 


Wind stirring. The across-front velccity patterns oenetratea 


j-- 


the mixed layer within the front and establish cross-fronital 
mixing. The ertire temperature field in the vicinity of the 
front undergoes an adjustment which inclines all the isct- 
hopmeedewawa=d tc the right in a rather uniform pattern. 

The mixed layér undergoes major adjustment over the 
G8-hcur integration, but the primary trends are a deepening 
SMepcen ssaqes cf the front, a leveling off on the right side 


agar aeeakaown cl the Shallow section within the frone due 


*o mixing and diffusion. By extrapolating the processes 
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Figure 31. Frent 2 Case II 48-Hour V-Velocity. Contour 
interval is 5 cm/s. 
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eccurrinag at 48 hours, it is presumed that the mixed layer 
depth would eventually even out with little or no remains of 
a frontal peak. AS in other cases, a marked deepsenirg in 


eoemPemed layer cccurs at the lef: boundary of «he front. 
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IX. MODEL RESULIS FOR FRONT 2 CA 


REL 


tte 


A. ANALYSIS AND DISCUSSICN 

Constant surface heating is combined with a ovoositive 
Wind stress. As the net Ekman transport wiil be from right 
to left, and with heating, we expect to see the formation of 
awartm plume c&£ water on the surface in the cool (left) 
ambient pcol. The surface heating should anhance «he inher- 
ent static stability and prevent the wind from compvleteiy 
Mixing the transrerted warm waters. This effect does hapren 
prior to the 12 heur poin= and becomes more pronounced over 
the entire integration (see Fig. 32). Downwelling and mix- 
ing cf the warned surface waters with the cooler undsrlying 


waters are cbhservesd at the leading adge of the progressing 


ct 
4 


Oe 


mie O 1) t iaweliing On he sight side of the ccoler wa 


Ii 
j-~ 


OP) 
a 
AD) 
th 


Wie eommeeplaces Water lost due £0 transport is cleariy 


+2 


t+ show 


(A 


Dy ccmparing successive temperature fields. The fro 
a tendency tc almost split at 30 m depth intc a Shallow sur- 
face front anda deeper frent. Perhaps a different combina- 


tion of forcing would make this happen. 
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T at hour 36 
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Figure 32. Front 2 Case Tit 36-Hour Temperature. Contour 
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The isctherme are much smoother than for the cases of 
wind stress alone. AS in Case I, the originally-horizcrtal 
field on «he isft sides becomes inclined, and the 


Originally-inclined field on the right side becomes level. 
The effect is net as dramatic as before. 

The mixed layer quickly responds to the applied heating. 
By 12 hours the mixed layer has re-formed in the far field 
at 4O m. The shallow section within the front is advected 
Witn the piume ct the front and mov2s into the cooler (ief*) 
eee fhe OLiginal a eeerer Section of the mixed Layer depth 
at the right tcundary of the ‘front has deepened to the mean 
Baee f2etad value frior to heur 12, Jel eingeen' at the rear 
as the shallow frontal peak is advected. This shallow "bub- 
ble" progress?2s to the left with the front causing shalicw- 


mixed 


iD 


ing in advance of the frort and deepening behind. Th 
layer responds with changes of up to 10 to 15 m ona time 
scale of semething less than 12 hours (compare Fig. 33 tc 
ba Gia 5) 26 

The along-frent velocity maintains a constant shape and 
GGnmetquzat2on over the 48-hour integration. The maxinrum 
speed decreases slowly over time from 100 cm/s at the start 


BOmolunecmyZ= DY hOUF 36. The horizontal gradient on the right 
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Figure 33. Frort 2 Case III 36-Hour Mixed Layer Depth. 
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side decreases rapidly but on thea left sida is preserved 
(Fig. 34) This is due to the advance of the plume and the 
resistance to mixing that the heating has established. 

At 24 and 36 hours, the across-front veiocity is neg- 


ative, with the maximum cf -10 cm/s occurring within the 


moving plume cf cverriding warm water (Fig. BS) es A wave 
oscillaticn in the v field is svident, and has the same 
direction and time scale as these of Cases I and II. Aeie 


MOUS, Veils pCSitive, which is probably a response to tns 
ageostrophic u-component created as the front undérgoes 
adjustment. At 48 hours the v-field has a small positive 


and Bay De inertial mction. 


Pee cemaanY ~~ CASE II 
The acplacaticn of heating i¢ able to generate and pre- 


serve a warm plume which ovérrides the cooler water in 


respense te net Ekman transport. Some dcownwelling is 
implied at che leading edge of the front. The lsotherms 
show a "kirk" and an ageostrophic u-velocity is there. Scme 


of the surface-warnmed water is mixed into the ambient pocl. 
Sn enes G2 gnc “uroscream"™ sids, upwelling occurs as cooler 


Tei eemracignc up to wecElace that lost by transpert. 


oF 





Depth (m) 


Figure 34. 


0.0 


25.0 


90.0 


100.0- 


125.0 


150.0 


175.0 


200.0 


225.0 


U at hour 36 


es | O 
‘ ‘ 
ee + 
ya : j 
“i . . ’ 5 
s ‘ " 
.ouiy : ; 
‘ w % ‘ ‘ 
Vit \ ; 
ee ne \ : 
ae ‘y , 
Sy ! 
Oo 1% 3 
1 f> 
i © t 
‘ : \ 
\ \ 
\ i 
\ ; 
\ _ ' 
\ \ 
\ ‘ ' 
\ ‘ ' 
\ me 
\ ve i 
\ ~~. ow 
\ 
\ 
\ 
\ 
\ 
\- 
\ 
\ 
\ 
\ 
\ 
\ 
‘\ 
x 
re 
ae 
iene ae a as 
25.0 35.0 45.0 Sanu) 65.0 


Y-distance (km) 


Beene 2) Caso iil 36>Houn U=Velocity. 


Con 
Sct ecme hemes 20 cm/s. 


5 
ct 
O 
f= 
It 


g5 





Depth (m) 


0.0 


25.0 


50.0 


75.0 


100.0 


125.0 


150.0 


175.0 


200.0 


225.0 


Vv amne@ur 24. 


rd : 
wit 
¢ 


ae me Sweets ete ee seas eer 
~~ 
cS 





SAX 


25.0 35.0 45.0 55.0 65.0 75.0 
Y-distance (km) 


Peegieeceemremertent 2 Case Lif 24-Hour V-Velocity. Contour 


Peecenval is 10 em/s- 


39 





The mixed iayer depth shallows to a constant depth outside 
or the frontal area under the uniform heating, and the shal- 
lower peak existent within the front is advected with the 
front in a resfonse to the mass transport. The across-front 
velocity generally shows a maximum at the surface within the 
front the direction of net Ekman transport, though reversals 
in dizrecticn cccur with time as a result of ageostrophic 


respenses. 
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ie ANALYSIS AND CISCUSSICN 

The cembined effects cf surface heating and a negative 
Wind stress act to push the upper 150 m of «he front inte 
the warmer water on the right side. Vertical mixing is pro- 
nounced as seen Ey the vertical orientation of the surfaced 
isotherms (Fig. 36). This result is expected as convective 
Mixing is ccupled with wind stirring. Diffusion and hori- 
zontal mixing cause the frontal horizontal temperature gra- 
dient to decrease continually over the 48-hour integration. 
Here *oo, the entire field develcps a downslope to the right 
aS in the preceding three cases. 

The mixed layer is initially shailowed +o around 40-m by 
the applied heating, and the shaliow section located within 
the front kreaks down and approaches the level of the cver- 


s expected based on 


weletexeGgmiayer (Fig. 37). This result 
previcus cas2s which combined the action of convective mix- 
fag and@waivd stirring. 

The u~velccity field haS an ageostrophic component at 


the left boundary aS a response to the adjustment. 
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Pigure 36. Frent 2. Case IV 36-Hour Temperature. Contour 
interval is 0.5 deg C. 
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Basically the u-velocity decreases its maximum speed and 


weakens its er 1. Zemea tL gradient over the 48-hour 


integration. The across-front v-velocity is in the direc- 


w 


ticn cf transport at 12 and 48 hours and extends frcem the 
surface, through the mixed layer depth, and down tc @ depth 
below 100 ao. OnegePadainyeccross-frontal mixing occurs. At 
24 heurs, the v-field is teversed over th2 entire upper 200 
Me + 36 hours, the v-efield exists only in the shallow sur- 
Pace Jayer, and, by 48 hours, the field returns to the pos:i- 
tive directicn. Some tyre of oscillation is evident in the 
Yemeescia (Zerecr to Figs. He, 39, CO and 41)< Ate houer 2. 
che field is fesitive everywhere. At howr 24, it is neg- 
ative everywhere. At Lour 36, it is negative in the shallow 
Gepths only. At 48 hours, 1% 1S G@gain positive everywhere. 
Minemrateern cectrred in the previous three cases, and it is 


scme type of oscillation in the 


rf 


strongly sugqc¢estive o 


V=3GQMEEOnenm= oCcth in direction and in depth-extent. 


B. Suman. -~ CASE IV 

The adtustrent of the frent under a negative wind stress 
and surface keating is affected by both convective mixing 
and wWand~sSstirring in the vertical. GzosS~frOrtale as x= ng Ss 


evident in the hcrizcntal. Micmracene 25 “oushed back! 2nve 


104 





Paemwacmer  Wacer Oo, the £ight by the transport of cocler 
water. The front diffuses over time. 

The nixed layer depth shallows cutside of the front due to 
the surface heating, and the shallow portion within the 
frent is ercded through mixing and approaches the depth of 
the far field. The u-velccity follows an uneventful pattern 
of diffusicn anda decrease in maximum speed as the front 
adjusts. The direction of the v-velocity correlates to the 
direction cf net Ekman transport only at 12and 48 hours. 
Between those times, an oscillatory pattern in directicn and 
depth is cbserved. This 1s the same pattern which was 
Socervedmin tke v-fieids cf Cases” 1, II, ana IIIf. The 
v+field indicates cross-frontal mixing as it extends over 4 


great depth. It is independent cf the mixed layer depth. 
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Figure 38. Front 2. Case IV 12-Hcur V-Velccity. Contour 
interval is 5 cm/s. 
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Figure 39. Frent 2 Case IV 24-Hcur V-Velocity. Contour 


interval is 5 cm/s. 
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Figure 40. fFrent 2. Case _ IV 36-Hcur V-Velccity. Contour 
interval is 5 cm/s. 
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Figure 41. Frent 2 Case IV 48-Hour V~Velocity. 
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MeeeeCONCLUSTOCNS AND REGCGUMENDATIONS 


en ee ee SE a = a eee ee eee ee Se 


A. GENERAL 

This thesis has examined the unsteady response of «we 
oceanic frents to the lccal atmospheric forcings of wind 
stress and surface heat flux. Ppeesihplvemry Of | maou 
enabled identification of btasic processes stimulated by the 
forcings, precessés which were not so clearly séén in Frent 
Die Front 2 was a mere realistic simulation of an oceanic 
DeoOnt,eaideenGtetULly, .t will serve as a stepping stone for 
EWGeE her Gnvestigations inte fronts and cesting of this fron- 
Sal model. 

There are a few very general properties that exist fer 
both Frents 1and 2, and which bear out our hyvootheses 
Sac Gime amchapter 2.8. 

1) Wand Gitection clays a primary role in frental 
adjustment. In the cases where warmer, less dénse water is 
transported <=c cccler, denser water, one Struceus one eons 
front is preserved. The mixed layer d¢pth is also preée- 
Saeved~eand Little cress~frontal mixing occurs. 

z) In the cases where cooler, denser water is trans- 
ported to warmer, less dense water, the horizontal tempera- 


ture gradient is diffused, the mixed layer depth within thes 
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front breaks dewr and deerens, and much cross-frontal nixing 
occurs. Cecnvective mixing cccurs along with wind stirring. 
2) A uniform constant surface heat flux shallows the 
mixed layer cutside of the frontal zone to a depth which is 
Peemerrlonab te wrbs amount cf heat flux. The combination of 
a +1 dyne/cme wind stress and a surface hsat flux of -0.004 


deg-cm/s mixes heat down toa depth of the initial mixed 


tt 


rh 


layer and leads te stratification there. Stieacir can yones 


ct 
- 


the isotherms at the surface did nox occur. For 1@ Cass 


) 


th 


in which tle heat flux is applied, turbulent mixing is 
restrained in the vertical, which agrees with the results of 
Niiler (1975) and Mellor and Durbin (1975). 

4) In ar¢cas of the frontal zone where ageostrophic 
velocities are created, vertical mixing is increased. 

5) Mixed layers deepen rapidly within che first 12 
hours under wind stress alcne as the fronts underge initial 
adjustment. Eeycnd 12 hours, adjustments are confined to an 
area around tke the basic depth attained during the initial 
adjustment period. Nigel cra oy Oe etlor and Durban(1975) 
and De Szceke (1980) all reported this phenomenon. 

6) A time-dependent pattern in the direction and 


@eprn-2nrtuence Of the across-fron= velocity exists. Tt ahs 








probably due tc inertial cscillations created by the atmcs- 
Bien c LObcinG. 

It was surprising to see <he downward bulge in «he mixed 
MeapeG at the cold-side boundary of both fronts in the wind 
stress only cases. It is unclear as to why this happens, 
though some speculative reasons were provided. 

In Front 1, the formation of cores of positive v-velicc- 
ity which are symmetric about the mixed layer depth are 
cra gua nd. Inertial oscillations may provide 2n 
explanaticn. 

The tilting down to tke right of the entire temperature 
field in Front 2 in all four cases is a major adjustment, 
and one which is independent of wind stress direction and 
heating. T= would appear that this tilt 1s a response +o 


eae re-orientation of the oceanic current fields. 


B. RECOMMENCATICNS 

There is Froth further investigation that can be done 
With the results of this thesis, as well 2s new work which 
2s suggested by these results. The results of these eight 
fedelesuns could te further quantified by using a finer con- 
Pommmemnecmyai fOr the fields, and by producing fields f£9r 


G@ead2ents, agecstrophic velocities, and the iike. The 
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inertial mcticns could be investigated to encompass nagni- 
tudes, directicns, depths end areas affected. t would be 
instructive tc ¢xamine {y,z) streamfunctions. The mcedel 
could be tested by conducting parametric sensitivity studies 
of the constarts used in the governing equations and mixed 
layer equaticnrs. 

New work might include using different combinations of 
Bhese atmeepheric forcings both in magnitude and direction. 
Inclusion cf a wind stress in the y-direction would cer- 
tainly prove to be of interest, as would a diurnal heating/ 
Gaolingeecyele. Cther atmcspheric forcings could be appiied. 
Evapcretion and precipitation might be incorporated. The 
equation of state could be expanded to include salinity, and 
an initial salinity field could be inserted and monitored 
over time. Cf course, a three-dimensional version of this 
model wceculd ke valuable. More simulations of observed 
fronts using data acquired in the field would previde fur- 
wcememcaam: into the applicability of the model. Data 
acquired by Jchannessen, e¢t al (1977) of the Maltese front 
would be a reccemmended starting foint. 

Computer central processor unit (CPU) time was extensive 


for each of these nodel runs. They each averagec over four 








Mewae OE GCkU tine +o Tun, EZoviding no stopping 
ing were required. SWcheamodeiecs far from 
operationally fer reali-tire assessment of ocean 
peev.ide, £6n Gxanols, input for acoustac foreca 


However, ftcr rasic research, this acde@i a 


AO, 
‘So 
a 
ey 
t 4 
W 


Sei acct onwa-c1phering the complexicvies of 


dynamics and thermodynamics. 


Se Vee =O 


sting acdels. 
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